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Abstract 



We model membrane proteins as anisotropic objects characterized by 



(N . 

^ . symmetric-traceless tensors and determine the coupling between these order- 

. parameters and membrane curvature. We consider the interactions 1) be- 

, tween transmembrane proteins that respect up-down (reflection) symmetry 

0> , of bilayer membranes and that have circular or non-circular cross-sectional 



areas in the tangent-plane of membranes, 2) between transmembrane pro- 
teins that break reflection symmetry and have circular or non-circular cross- 



O . sectional areas, and 3) between non-transmembrane proteins. Using a 



field theoretic approach, we find non-entropic interactions between 

reflection-symmetry-breaking transmembrane proteins with circular cross- 
sectional area and entropic interactions between transmembrane pro- 
teins with circular cross-section that do not break up-down symmetry in 
agreement with previous calculations. We also find anisotropic interac- 
tions between reflection-symmetry-conserving transmembrane proteins with 
non-circular cross-section, anisotropic l/i?^ interactions between reflection- 
symmetry-breaking transmembrane proteins with non-circular cross-section, 
and non-entropic l/R"^ many-particle interactions among non-transmembrane 
proteins. For large R, these interactions are considerably larger than Van der 
Waals interactions or screened electrostatic interactions and might provide 

1 



the dominant force inducing aggregation of the membrane proteins. 
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I. INTRODUCTION 



Recently, the structure and properties of model and biological membranes have been 
studied extensively. Biological membranes play a central role in both the structure and 
function of cells. Biomembranes divide living tissue into different compartments or cells and 
act as cell boundaries. They determine the nature of all communication between the inside 
and the outside of cells. This communication can take place via the actual passage of ions or 
molecules between two compartments or via conformational changes induced in membrane 
components. Model bilayer lipid membranes in aqueous environments exhibit many of the 
attributes of the biological membranes. For example, these membranes can form vesicles 
or more complex structures that divide space into separate compartments, which like cells 
can fuse or divide. However, there are many properties of biomembranes that cannot be 
mimicked by lipid bilayers. Energy-driven transport of ions across membranes and receptor- 
mediated events are only a few of the myriad of membrane-associated functions that lipid 
bilayers are incapable of performing on their own. Such processes are mediated by proteins 
that are attached to or dissolved in biological membranes Thus in order to make 

lipid bilayers more realistic models of biomembranes, it is necessary to introduce into them 
membrane inclusions such as proteins. 

Membrane proteins are classified as integral proteins or peripheral proteins according 
to how tightly they are associated with membranes. Integral proteins are so tightly bound 
to membrane lipids by hydrophobic forces that they can be freed only under denaturing 
conditions. Peripheral proteins associate with a membrane by binding at the membrane 
surface; they can be non-destructively dissociated from the membrane by relatively mild 
procedures. Some integral proteins, known as transmembrane proteins, span the membrane, 
whereas others are attached to a specific surface of a membrane. For brevity, we refer the 
latter proteins as non-transmembrane proteins to distinguish these from transmembrane 
proteins. However, no proteins are known to be completely buried in a membrane 0. 

Interactions between membrane proteins is expected to be controlled by lipid affinity. 
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direct interactions such as electrostatic and Van der Waals interactions, and indirect inter- 
actions mediated by the membrane 0. The latter interaction arises from thermally-driven 
undulations of the membrane and is analogous to the Casimir force between conducting 
plates. Since the degree to which fluctuations of a membrane are restricted depends on dis- 
tance between membrane proteins, its free energy also depends on this distance, decreasing 
with decreasing separation. This implies an attractive force, which leads to a tendency for 
membrane proteins to aggregate. This indirect interaction between membrane inclusions 
was first calculated by Goulian, Bruinsma, and Pincus 0. Before presenting the interaction 
models introduced in this paper in detail, in Sec. II we will review its results, indicating where 
they differ from and extend those of Goulian et al. In Sec. Ill, we introduce a phenomeno- 
logical model (Model I) for protein interactions. In this model, proteins are characterized by 
symmetric-traceless tensors depending on the cross-section shapes of proteins on the mem- 
brane, and interactions are described by symmetry-allowed couplings between these tensor 
order-parameters and the curvature tensor of the fluctuating membrane. In Sec. IV, we intro- 
duce another phenomenological model (Model II), in which there is an interaction between 
a membrane protein and membrane lipids at its perimeter. At the perimeter, lipids tend to 
align with the direction of the protein at a certain angle depending on whether or not pro- 
teins break the up-down symmetry of the bilayer membrane. The interaction is described by 
the fluctuation of the normal vector of the membrane around this preferred direction. The 
interaction between non-transmembrane proteins is described in Sec.V. Non-transmembrane 
proteins are proteins with preferred center-of-mass positions not at the center of the bilayer 
membrane. We expand the potential energy in terms of the deviation from this preferred 
position and include other couplings considering the symmetry to calculate the interaction 
between non-transmembrane proteins. Finally, a discussion is given in Sec. VI. 



II. REVIEW AND SUMMARY 
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A. Review of the previous work 



The interaction between membrane inclusions with circular cross-section on the mem- 
brane has been calculated by Goulian et al. . They use the Helfrich-Canham Hamiltonian 

m = \j d\^{nH^ + RK), (2.1) 

to describe fluctuations of the membrane, where k and R are the bending and the Gaussian 
rigidities and H and K are the mean and the Gaussian curvatures, respectively. The surface 
tension is not taken into account since it effectively vanishes []TU|,|rT|. Within inclusions 
with circular cross- sect ion, the constants k and R are assumed to differ from those of the 
surrounding membranes. For some inclusions, such as proteins, the circular regions are 
assumed to be rigid with n = —R = oo. This case is the strong-coupling regime. On the 
other hand, for regions with excess concentrations of lipids, /t and R within circular regions 
are assumed to have values close to those of the surrounding membrane. This case is the 
perturbative regime. In both the strong-coupling and perturbative regimes, Ref. |^ flnds 
that there is an entropic l/R^ interaction, which is proportional to the temperature and 
to the square of the area of the circular region. Also, at low temperatures, there is a non- 
entropic 1/ interaction between proteins which varies with the square of angle of contact 



between membrane and proteins. Recently, Golestanian, Goulian, and Kardar |T2[ extended 
the calculation in Ref. |0] to the interaction between two rods embedded in a fluctuating 
membrane. They flnd an anisotropic l/R^ interaction between two rods. 



B. Summary of the present work 

In this paper, we introduce three models for the interaction between membrane inclu- 
sions such as transmembrane proteins and non-transmembrane proteins. First, we present 
Model I. Here, proteins are characterized by symmetric-traceless tensor order-parameters, 
and the coupling between these order-parameters and membrane curvature is determined 
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by symmetry and power counting. If proteins don't respect up-down symmetry, we allow 
for couplings that break up-down symmetry. Otherwise, we require up-down symmetry. 
For proteins with circular cross-sectional area on the membrane and preserving up-down 
symmetry, we find the leading distance-dependent free energy. 



^ = A 2 2 04 / / ^ ^ Mx)SK{y) + 6K{y)6R{x)] , (2.2) 

where the coupling constants are denoted in terms of the bending and the Gaussian rigidity 
differences, 5k{x) and 5r{x), between proteins and surrounding membrane. This corre- 
sponds to the perturbative regime of Ref. 0, whose results we reproduce. For proteins with 
circular cross-sectional area, up-down-symmetry-breaking couplings do not affect the lead- 
ing contribution to the free energy, and the leading distance-dependent free energy remains 



identical to Eq. (2^). We also calculate interactions between proteins with non-circular 



cross-sectional area. When up-down symmetry is conserved, we find an interaction energy 

^ = -^g^ 64^2^2 j:g4 ((^4^4 + q2Ql){q4QA + q2Ql + d2Ql) cos 4(^1 + 62) 

+2(^4^4 + q2Ql + d2Ql)d2Ql cos^ 2^1 cos^ 262 + {q2Ql + d2Ql)q2Qt) , (2.3) 



where A is the cross-sectional area of proteins, Q2 and Q4 are magnitudes of 2nd-rank and 
4th-rank tensor order parameters measuring orientational anisotropy, respectively, and 9i 
are the angles of the directions of proteins measured with respect to the separation vec- 
tor between proteins (See fig. [l|). This anisotropic l/R^ interaction contains anisotropic 
cos^ 29i cos^ 26*2 interaction in addition to anisotropic cos 4(0i + 62) interaction also found in 



the recent independent work by Golestanian et al. |T^. However, when up-down symmetry 
is broken, there is an anisotropic 1/B? interaction: 

Thus the leading term in the free energy falls off with separation as rather than 1/ R^. 

Next, we introduce Model II in which we impose a certain boundary condition at the 
perimeter of proteins with the circular cross- sectional area. For proteins with up-down 
symmetry, we find 
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T = -kuT 



6^2 



(2.5) 



This free energy looks similar to that of Ref. 0] in the strong-coupling regime. When proteins 



break up-down symmetry, we find 



6^2 



(2.6) 



+ 



where is the contact angle between the direction of i-ih. protein and the unit normal of 



Finally, we introduce a height-displacement model in which protein positions normal to 
the membrane can vary. In this model, we find three- and four-body interactions in addition 
to a two-body interaction. These three- and four-body interactions also fall off as and 
are the same order of magnitude as two-body interaction. 

Consequently, by introducing three models to describe the interaction between membrane 
inclusions such as proteins, we recover all the results in Ref. 0. Furthermore, we obtain 
anisotropic interactions between proteins with non-circular cross- sectional area. Also, we 
extend the calculation to the up-down symmetry breaking proteins with non-circular cross- 
section and find anisotropic interaction between them. Moreover, using a height- 
displacement model, we find three- and four-body interaction in addition to two-body 
interaction. 

III. MODEL I 

For a fluid membrane free of membrane proteins, the energy of membrane conformations 
can be described by the Helfrich-Canham Hamiltonian [H,^, 



expressed in terms of the local mean and Gaussian curvatures. We will work at length scales 
large compared with the membrane thickness but small compared with the membrane's 
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the membrane (See fig. §). In the limit T — > 0, this free energy becomes the result in Ref. 
for the low-temperature regime. 




(3.1) 



persistence length. Thus, we can parameterize the membrane in the Monge gauge R = 
(x, /i(x)) where x = {ui,U2)- In terms of R and the unit normal vector of the membrane 
N, the metric tensor gap is given by S^R • S^jR and the curvature tensor Kap is given by 
N • DaDjjR, where Da is the covariant derivative along Ua direction on the membrane. In 
the Monge gauge, 

H^^g^^Kaf,^V^h + 0{h'), (3.2) 

K = det g'^^Kp^ = V^hV'^h - dadphdadph + (3.3) 

When the topology of the membrane is fixed, the Gaussian curvature term can be dropped, 
and the leading term in TYq in an expansion in derivatives of h is 

Ho^^kJ (fxW^hW^h. (3.4) 

Now let us consider the coupling between membrane proteins and membranes. 



A. Proteins with circular cross-section 

Membrane proteins can have arbitrary shapes; as a result their tangent-plane cross- 
sections can be any shape. Now we will compute the undulation mediated force between 
proteins separated by a distance larger than the size of proteins. For simplicity, let us 
first consider membrane proteins that have a circular cross-sectional area on the membrane. 
These proteins may be described by a scalar density, p, which may be interpreted as the 
distribution function of proteins describing the positions of proteins and the configurations 
of protein's amino acid sequence 

P(x) = I]^/i(x-Xi), (3.5) 

where x = {u^, u^) is a point on the membrane and the sum is over all proteins. The specific 
form of /j(x — Xj) depends on the specific conformation of i-th protein at the position x,. It 
vanishes outside the protein cross section: 
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f{x), |x| < ap, ^^^^ 



0, |x| > ttp. 

where Up — ^JA/tt is the radius of the protein where A is its cross-sectional area. We assume 
all proteins are identical so that they are all described by the the same function /(|x — Xj|). 
For membrane proteins, ap is of order lO^A. If we model the protein as a uniform cylinder, 
the distribution function of protein will be f{x) = 1 inside the projected area and f{x) = 
outside. In general, proteins have non-uniform folding of the amino acid chain, and f{x) 
will have small deviations from unity inside the circular cross-sectional region D. In this 
case, we use 

/ (fxfix) = A (3.7) 
Jd 

as the definition of A. 

When proteins do not break up-down symmetry, the relevant coupling between p and 
the height fluctuation fleld of the membrane is 

Hint d^x[ap{^)KlKl + ip{^)KlK'^ 

-\T.I dFx[af{\^-^\)KlKl + ^f{\^-^\)KlK% (3.8) 

where Dj, i = 1, 2, • • • denote circular regions occupied by membrane proteins. The coupling 
constants a and 7 describe couplings between the density inside protein's cross section and 
the curvature of a membrane. Thus these can be related to the bending and Gaussian 
rigidities: 

af{x) — Sk{x) + Sr{x) , "yf{x) — —Sr{x) , (3.9) 

where Sk and SR can be interpreted as the changes in the bending and the Gaussian rigidities 
due to the existence of proteins on the membrane. In the Monge gauge, to lowest order in h 

Kl = -V^/i , K^K^ = dad,hdM, (3.10) 

and the relevant coupling becomes 



Hint = lJ d^x[ap{^){V^hf + ^p{^)dad,hdad,h]. 
The free energy is given by 



(3.11) 



-^P / d''x[ap{:K){V^hf + -ip{^)dad,hdadbh]. 



(3.12) 



We can use the cumulant expansion to calculate this form of the free energy. We write 

^ ^exp[-^/5 j d?x[ap{i^){V^hf + ^ p{i^)dad,hdad,h]]^ ^ (3.13) 
where ()o denotes the ensemble average over the fluid membrane Hamiltonian only and 

^-m = j [p;,] exp[-^/5fi: j d^x{V^hf]. (3.14) 
The cumulant expansion gives 



(e% = (l + V^+-V^^ + ---)o 



exp 



{V)o + ^[(y\~{Vf^ + 0{V^) 



(3.15) 



Plugging Eq. ( |3.13| ) into the cumulant expansion Eq. ( ^.15| ) and keeping terms up to order 
/i^, we find the free energy 

-/5(-F-^o~ ' ^ 



2-/3 j d'x[ap{^)iy'hY + -ip{-K)dad,hdad^h]^ ^ 



~\ (^-\l3 j d''x[ap{^){y^hf + ^p{^)dad,hdad,h]^^ . (3.16) 



This can be expanded in terms of the height correlation function Ghh{^ — y) and its deriva- 
tives. The height correlation function in the real space is 



Ghh{y^-y) 



H^)Ky))o 

(27r)2 (27r)2 
d'^p e*P (^-y) 1 



(27r)2 /3kp^ 
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{Hp)hiq))o 



IGtt/^k 



(3.17) 



where we used (^(p)^(q))o = (2vr)^5(p + q)//3Kp^ in momentum space and x — y = R. 
Then, taking derivatives we find 

{dadbh{x.)didjh{y))o = dadbdidjGhhi^ - y) 



AtxI3kR^ 

— 2{RaRbSij + RaRiSbj + RaRjSbi + RbRiSaj + RbRjSai + RiRjSab) 



-SRaRbRiRj 
1 



(3.18) 



where R = R/i?, R = |R|. We proceed to calculate the terms in Eq. ( ^.Ifil ). We are 
only interested in terms that depend on the distance between membrane proteins. We can, 
therefore, drop the first and the last terms in the RHS of Eq. (|3.16| ) since they do not depend 
on distance: 



J (fxp{-x)dadbhdadbh^^ = J (fxp{x)V^Ghh{^-y) 



iy=x 



(27r)2/5 



(Pp 1 



N 



(27r)2/5« P 



N 



;a,A)^ 



and similarly 



[apAf 



(3.19) 



(3.20) 



4/?K 

where is the number of proteins. In Eq. (p.l9|) we introduced the cut-off for the height 
fluctuation fleld A ~ l/a^ where ap is the radius of the protein. The second term gives 
contribution to the distance-dependent free energy: 



^ /rfWyp(x)p(y) 4a79,9,V^G,,(x-y)9„AV^G,,(x-y) 



+2-f^dadbdidjGhh{^ - y)dadbdidjGhh{y^ - y) 



(3.21) 



^ / d'^xd'^yp{x.)p{y)- 



4a-fTabii{R)Tabn{Il) + 2-f^Tabij(R)Tabij{R) 



{An/SKR^y 

where we kept only the leading distance-dependent terms and R = x — y. Thus, the leading 
distance-dependent free energy is given by 
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-(3r=\jd?xj rf^yp(x)t/(|x - y|)p(y), (3.22) 

where 

For two proteins separated by a distance i?, the leading R dependence is 

- = M L ^'^ L ["-' ^ 

Relating the couplings a and 7 with the variations of the bending rigidity and the Gaussian 
rigidity as in Eq. (|3.9|), we recover the result of Goulian et al. M 



^ = t¥t^a I / [^<^)^'<y) + 5<y)5-<x)] . (3.25) 
47r^K^/i* JDi JD2 

If membrane proteins break up-down bilayer symmetry, there is another possible relevant 
coupling, 

= \\j d'xpi^)K: = IxY: rfV(|x-x,|)ir:. (3.26) 

However, this term does not contribute to protein-protein interactions since the distance- 
dependent contribution vanishes as follows, 

ci2a;rf2^p(x)p(y)(V2/i(x)V2/i(y))o = J dVypi^)p{y)\/''Ghh{^ - y) (3.27) 

= 2 r r d'yf{x)f{y)6{x~y) = 0. 

(dK JDi JD2 



B. Proteins with non-circular cross-sections 

So far we have, for simplicity, considered protein-protein interactions when proteins have 
circular cross section. However, in general proteins have asymmetric conformations giving 
rise to non-circular foot prints on the membrane surface. They can then be characterized 
by symmetric-traceless tensor order-parameters such as Q"*? Q"''"^'^, and so on: 



12 



Q'^'i^) = E - X,), Q'^'-^'i^) = E -^QT'fi^ - ^m), (3-28) 

where denotes the position of protein and and Q'^^'^'^ are the symmetric-traceless 
tensors constructed from the characteristic direction vector of /^-th protein on the membrane. 

When up-down symmetry is not broken, the relevent couphng between inclusions and 
curvature is 

Wint d^xS'''^(x)Kaj,K,a, (3.29) 

where 

^abcd^^) ^ 54Q"''"'(X) + g2Q"''(x)Q"'^(x) + d2g"^(x)g^''(x). (3.30) 

The couphng constants q^, q2 and ^2 describe couphngs between the protein tensor order 
parameters and the curvature of a membrane. Results for membranes with circular cross- 
sections can be obtained by choosing 

Sf^ = aS^^S"^ + ^5'^5^'^, (3.31) 

rather than insisting the order parameters be symmetric and traceless. With this coupling, 
we proceed as before using the cumulant expansion. For two proteins separated by a distance 
vector R from one to the other, the free energy becomes 

-(5T='^ j d^j A^"'''(x)5'^"(y)Rafe/i(x)aead/i(x)a,a,/i(y)afeai/i(y))o (3.32) 



d^x j d^yS-''^{^)S'^^\y)dad,did^GnH{^-y)dad,d,d^^^^^ 
In terms of the tensor TaUj introduced before, the final form for the free energy writes as 

-(3j^^-^j dFx j d^yS^"^{^)Uai,ai,m{\^ - y\)s''''\y), (3.33) 



where 



jT /I , IN {TaUj{fCjTcdkl{^) + Tai,kl{fi)Tcdij{B.)) 

Uabcd,ijkl{\^ - y|J = g^2^2|x_y|4 " ^'^''^^■l 
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For two identical proteins separated by R = Xi — X2, the leading distance-dependent free 
energy is found to be 

A2cabcdf \cijkU \ 

^ ^ -^^^ 647rV2^4 [Tabij{^)Tcdki{^) + TMR)ndijm. (3.35) 
Now the free energy is anisotropic, depending on the direction of the separation vector R and 
the orientation of proteins described by S"''"^{xi). Q"-^"^ is a 4th-rank symmetric-traceless 
tensor, which can be expressed as 

6 

^^^^ab^cd ^ ^ac^hd ^ rt''")) , (3.36) 

and (5"* is a 2nd-rank symmetric-traceless tensor; 

Q'^' = Q2 {e'^e' - ]^6'''^ , (3.37) 

where = cos 6 and = sin 9 characterize the direction of protein with 6 measured with 
respect to the separation vector R and Q2 and (^4 are magnitudes of 2-fold and 4-fold 
anisotropy, respectively. Then, the free energy becomes 

^ = -fei?^ g4^2^2^4 ((^4^4 + q2Ql){q4Q4 + ^2^2 + ^^2^2) COs4(^i + 62) 

+2(?4<54 + g2Q2 + d2Ql)d2Ql cos^ 2^1 cos^ 2^2 + (?2<52 + d2Ql)q2Qt) , (3.38) 
Again, for proteins breaking up-down symmetry, we have the additional relevant coupling 

Hint = lj d^xS^\^)Kab. (3.39) 

In contrast to the case of circular cross section, this coupling leads to a qualitative change 
in the protein-protein interaction. Proceeding as above, we find the leading distance depen- 
dence of protein inteaction is 

-Pr=^jd^xj d'yS''\^)S^^{y){dad,h{^)d,d,h{y))o 
-^Jd'xJ d^yS''\^)S'\y)dad,didjGHh{^-y) 
-\jd'xj A5"^(x)t/„,,,(|x - y|)5*^(y), (3.40) 
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where 



Uab,ij{\^-y\) 



(3.41) 



27rK|x — y 



This interaction is also anisotropic, depending on R and S''^^(xi). For spherical cross section, 
since S""^ ~ 6ab and d'^^d^^Tabij = 0, the contribution to the interaction vanishes as before. 
For ellipsoidal cross section, S""^ = diQ"'^ where di is the coupling constant and the free 
energy becomes 



The minimum energy configurations are at 6i + 62 = 0,7^. 

Consequently, by introducing symmetric-traceless tensors as the order-parameters for 
anisotropic proteins and by determining the relevant couplings by symmetry, we were able 
to rederive the results for the circular cross section by Goulian et al. [0. Furthermore, we 
obtained anisotropic interactions between proteins which have the non-circular cross section. 
This anisotropic interaction has the leading distance dependence l/R^ and l/i?^ depending 
on up-down symmetry breaking. 



In the previous section, we introduced a coupling between membrane proteins and the 
height fluctuation fleld of the membrane by considering symmetry and power counting. 



distribution function of proteins, the physical implication of this coupling can be that the 
bending and Gaussian rigidities inside the protein cross section differ slightly from those of 
the surrounding membrane. Thus this coupling can be thought of as perturbative. However, 
if proteins are infinitely rigid with n = —R = 00 inside the protein cross section, perturba- 
tion theory fails. In this case, we can derive the protein-protein interaction by considering 
the phenomenological interaction between membrane proteins and membrane lipids at the 
perimeter of the proteins. 




(3.42) 



IV. MODEL II 



Since the order parameter for proteins in the coupling Eq. ( p.8|) can be interpreted as the 
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First, let us consider proteins which have circular cross sections and do not break up- 
down symmetry. These proteins can be modelled as inversion-symmetric three-dimensional 
ellipsoids of revolution (or cylinder) with a major axis pointing along a unit vector m in 
three-dimensions. Their orientational order can be characterized by the symmetric-traceless 
tensor Qij = {rriimj — Sij/3). We assume the axis m prefers to align along the membrane 
normal N. A simple interaction favoring this alignment is 

^mt = -^«^/ -^g^^iv.iv,. (4.1) 

In the Mongc gauge, 

Q'^NiNj = (m ■ N)2 - 1 = -S'^'ima - Na) {nib - N^) + constant, (4.2) 

where a, b run over 1,2 only and Na — —dah to lowest order in h. Now we can Taylor-expand 
the unit normal of the membrane at the perimeter from the center of the protein to lowest 
non-trivial order in Up 

Na{r)\c, = iV„(R^ + aph{l - (b • m^)Y^) = 7V„(R^) + a^b • V7V„ + • • • , (4.3) 

where b is the unit vector from the center of the protein to its perimeter and A^a(R^) is the 
average of Na,{r) along the perimeter c^. Dropping the constant term, the coupling becomes 



r dl 

'Hint = a II / t; [ma(R^,) - Na(R^) - aph^d^Na] 

„ Jca iTldn 



p 



m„(R^)m„(R^) + -^d,Nad,Na 



where ma(R^) = ma(R^) - NaiRij,). The free energy is 



(4.4) 



e-^^ = J [P/i][Pm„(R^)]e-^^°-^^'"* (4.5) 
= (constant) J [Dh] exp[~(3K J d'^xiy^hf - ^l3aAY.dbNa{B.^)d,Na{B.^)], 

where Hq — J (PxiV'^hY and the integration over rha is trivial and gives constant contri- 
bution. Since the coupling has a quadratic form, we can evaluate this using the Hubbard- 
Stratonovich transformation. Although it is nothing more than completing the square, we 
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will find this technique to be very useful. By introducing the auxiliary fields W(/i) and 
defining V(/i) as 

V,{i^) = d,N,{i^) , V^ii^) = V2d,NM = V2d2m{ij.) , VM=d2N2{i^) , (4.6) 
we have 

e-P^ = f [m][DW(/x)]exp[-/37^o - ^Y.^i^') ' W(/.) + ^ ^ W(/x) ■ V(/x)] 

= e-^^^ ( [VW{^^)] eM-^ E W(/i) ■ W(/i)] (exp[^ ^ W(/i) ■ V(/i)])o, (4.7) 
where F = j3aA/'n. Using the cumulant expansion again, to the lowest order we obtain 
e-^("«) = / [I^W(^)] exp[-( V W(/.) ■ V(/.)W(z.) ■ V(z.))o] 

= /[PW(/x)]exp[-i5:H^„(^)(K(^)H(^))oW^6(^)]. (4.8) 
For two proteins separated by a distance i?, we find 

e-^("«) = /[PW(l)][I?W(2)]exp[-l ^ PFa(/i)(K(/i)n(z^))oW^.(^)] 
= (det(K(/x)H(^^))o)"'/' 

= (constant) fl - 3(-|-)2') , (4.9) 



where A ~ l/a^ is a cut-off for the height fluctuation, and is the radius of the protein 
introduced in Sec. II. A. Thus, the free energy has an i?-dependence as 

96 ^J? 

in accord with the previous calculation by Goulian et al. [1^ . In the above equation, we used 



a cut-off for the height fluctuation, A = 2/ap [jT3 



For proteins that break the up-down symmetry, the unit normal of the membrane at the 
perimeter of the protein is not forced to be parallel to the direction of the protein. Instead, 
the unit normal is forced to have a fixed angle with the direction of the /i-th protein. 
Thus the preferred unit normal at the perimeter is 
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No(r)|c, 



m + a^b 



(4.11) 



The coupling between the protein and the membrane at the perimeter of the protein is 

dl 

We Taylor expand to find to lowest order 



(4.12) 



r dl 



(4.13) 



(4.14) 



The free energy is now 

e-^^ = I [Vh][V7ha{i^)] exp[-^(3K J d^x{V^hf 

= / [m] exp [ d-'xiV^hf - ^ J2[9bNadbNa - '^dM 

J \_ Z J ZTT ^ ttp 

Re-defining V(//) as 

V,{i^) ^ d,N,{i^) - ^ , V^ii^) ^ d,NM ^ d2N,{i^) , VM = d2N2{i^) - ^ , (4.15) 

Clp dp 

the free energy becomes 



jm 



exp 



\f^Kjd'x{v'hf-^Y.^{^,).Y{^,) 

= |[m][PW(//)]exp -l/5«:|c^^x(V^/i)^-^EW(/i)-W(//)+iEW(//)-V(//) 

= / [PW(/.)] exp[-i; E W(/i) • W(/.)](exp[z E W(/.) • V(//)])o, (4.16) 
where F = PaA/n. For two proteins separated by a distance we find 
e-/^^ = / [DW(l)]pW(2)] exp[^ ^ W(/.) • (V(/.))o 

-\ E ((W(/.)-V(/.)W(i.)-V(i.))o-(W(/.)-V(/.))o(W(i.)-V(i.))o)] 

\WMW-iAW 



(detM)-V2e-|^^-'^, 



(4.17) 



where 



W = (W^i(l), 1^2(1), W^3(l), W^i(2), W2{2), 1^3(2)), 
^=(—,0,^,^,0,^), 



M 



A2 



327r/?K 



3 1 -A -A 

2 -2A 

1 3 -A 3A 
-A -A 3 1 

-2A 2 
-A 3A 1 3 



(4.18) 
(4.19) 



(4.20) 



with A = 8/ i?^ and A = 2/ap is the membrane cutoff. Thus the i?-dependence of the free 
energy becomes 

1 



2^ 1 IGtt/Sk^^ ^ 1 



I ( D"? A 2 ) 



2 A2 



0(af + a^)/< 



2'i?2A2 

Our final form for the free energy is 



a2A6 



(4.21) 



J" 



vr 



(4.22) 



This gives the previous result Eq. ( [4.10| ) for = which corresponds to the strong-coupling 
regime in Ref. 0. In the limit T ^ 0, this gives the result for the low temperature regime in 
Ref. 0. Thus, in this phenomenological model, we obtain the general interaction between 
the up-down symmetry breaking proteins at finite temperature T. 

This calculation, which focuses on the change in free energy brought about by the ad- 
dition of inclusions, does not show explicitly how these inclusions modify the shape of the 
membrane at large distances from the inclusions. Careful treatment of the minimum energy 
configuration of h, about which we calculated Gaussian fluctuations, yields the same large 
distance distortion {h ~ cos n6') as calculated by Goulian et. al. We believe this result to be 
true for a free membrane with no imposed boundary conditions. If the membrane is forced 
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to be flat by an aligning field, these long-range forces will become short-range and surface 
tension will have a similar effect. It is not so clear what will happen on a vesicle of spherical 
topology with no Laplace pressure. This question is currently under investigation. 



V. HEIGHT-DISPLACEMENT MODEL 

Non-transmembrane proteins are exposed to a specific surface of a membrane. Thus, 
they have preferred center-of-mass positions not at the center of the bilayer (See fig. ^). We 
consider the interaction between these proteins by introducing the potential energy '^{( — h) 
where ( is the position of the protein and h is the membrane height fiuctuation field. For 
integral proteins, — /i) has a minimum at the non- vanishing value of ( — h = tq. We 
can expand — h) in terms of the deviation from this preferred value tq 

vl/(C -h) = ^{ro) + ^((C -h)- ro)2^|(c-.)=.o + • • • , (5.1) 

and if proteins are tightly bound, ^^\((;-h)=ro ^ 1- Considering the symmetry, we introduce 
the couplings 

^int = ^ E / d\[h{in + Vhf + k2{C- h)V^h 

+A;3(C - h)mambdadbh + ^(C - h)], (5.2) 

where m is the preferred direction of the protein. By minimizing over ma, we find 

ma = -dah + ^(C - h)d,hdad,h + Oih"). (5.3) 
ki 

Substituting this result into the coupling, we obtain in lowest order 

^int = ^ E / d'^AHC - h)V^h + k^iC - h)dahdthdadth 



+vl/(ro) + -((C -h)- ro)2^|(c-.)=n,]. (5.4) 



Minimizing over ((^ — h) gives the preferred position of the protein as 



(-h = ro-k2 (^\i(-h)=ro] V^/i + Oih'). (5.5) 
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Thus we obtain the couphng 

Hint = ^ E d^x[hroV^h + hrodahd,hdad,h - hi |^^|(^_,)^,^ j [y^hf]. (5.6) 

The first and the last terms look similar to the ones in the phenomenological model. However, 
the non-linear second term is allowed because the up-down symmetry is broken by the 
preferred position of the protein. In the low temperature limit, we assume lipids are so 
tightly bound to the proteins that (^■|(c-/i)=r'o) is much bigger than k2 and we drop the 
last term in Eq. (|5.6|) . In this limit, the Hamiltonian becomes 

n = ]-K I (fx{V^hf + J V / d^x[k2roV^h + ksrodahdbhdadbh]. (5.7) 
2 J 2 ^ Jd^ 

Now by minimizing this Hamiltonian over h, we obtain the low temperature limit for the 
protein interactions. From the minimum condition, 

^ = = KV'h + Y: k2r,VH{r - r,) 
6h y 

+ E hroC^daSir - r ^)di,hdadi,h + dahdbhdadb5{r - r^)), (5.8) 

we obtain the equilibrium height for the membrane 
h{r) = Y^G{v-v^)-^ I Sr'G{r - r') ^ 9„G(r' - r^)9,G(r' - r,)a,a,G(r' - r,). 

(5.9) 

Thus the interaction becomes 

n = -h fd'x (j2 ^^(r - r^) + ^ ^ ^.^(r - r,)dbG{r - r,)dadbG{r - rx)] 

E f-^) - T,)dbG{v^ - T,)d,dbG{r^ - r.). (5.10) 

For two proteins separated by a distance i?, we obtain the leading distance dependence of 
the interaction 
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n = 18hrodaG{R)dbG{R)dadbG{R) 



-IShro 



271K / 



fk2roYl8k3ro , . 

= [2^)^- 

Above we used G{R) = — (A;2?"o/4vrK) InR^. We can interprete the parameters in Eq. (^.111) 
as the area of proteins and the contact angle between proteins and hpids. 

When there are several proteins, from Eq. ( ^.lOj ), we find that three- and four-body 
interactions exist in addition to two-body interaction. For three proteins separated by R^i, 
which is the vector from the u-th protein to the /i-th protein, we find three-body interaction 
to be 

where J2' means all fi, u, a are different, in addition to two-body interaction between each 
pair of proteins given by Eq. (|5.11|) . Similarly, we find four-body interaction to be 

where Yl' means all /i, u, A, a are different. Note that these three- and four-body interactions 
are also interaction which is the same order as two-body interaction. 



VI. DISCUSSION 

We model biological membrane as a continuous bilayer of lipid molecules in which various 
membrane inclusions such as proteins are embedded. Such model membranes with inclu- 
sions also have potetial applications for target drug delivery, nano-scale pumps, functionized 
interfaces, and chemical reactors. In this paper, we study how the membrane contributes 
to the interactions between inclusions. Also, it is interesting to understand how inclusions 
affect the properties such as rigidity or shape of model membranes. 

The interaction between membrane inclusions such as proteins with circular cross- 
sectional area was first calculated by Goulian et al. Using three models, which we refer 
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to as Model I, Model II, and a height-displacement models, we recover all the results by 
Goulian et al. The interaction in Eq. (2.2) is a temperature-dependent interaction between 
two circular inclusions that falls off with distance as 1/ R^. Assuming two inclusions are iden- 
tical, the force will be attractive if 5k and 5R have the opposite sign and repulsive otherwise. 
For inclusions with up-down symmetry, the interaction is attractive and falls off as 
again. The magnitude is set by ksT and is independent of the rigidity k. When inclusions 
break up-down symmetry, in addition to the attractive interaction set by ksT, we find a 
repulsive interaction proportional to the square of the contact angle in Eq. (2.6). Thus, for 
up-down asymmetric inclusions, there are competing attractive and repulsive interactions 
and we might have an interesting transition between aggregation and mixing of inclusions 
when ksT ~ ko^. Both Model I (for soft inclusions) and Model II (for hard inclusions) 
predict potentials that fall off with distance as R~*. This interaction is attractive for hard 
inclusions. For soft inclusions the sign of the interaction depends on the relative sign of 5k 
and 5k and is attractive if they have opposite signs. Reasonable models predict 5k5k < 
so that the prediction of models I and II can be viewed as being consistent. 

Furthermore, we calculate the interaction between proteins with the non-circular cross- 
sectional area and find anisotropic and interactions depending on whether up- 
down symmetry is broken or not. In Eq. (2.3), we find the interaction between proteins 
with non-circular cross-sectional area when up-down symmetry is conserved, and the free 
energy again falls off as 1/ but the magnitude depends on the orientations of inclusions. 
When only 4-fold anisotropy is nonvanishing, the angular dependence of the interaction is 
of the form cos4(^i -|- 62) ■ The interaction depends on the relative orientations of inclusions 
to the separation vector from one inclusion and the other. The interaction is attractive if 
1^1 + 6*21 < 7r/8 and repulsive if 7r/8 < 16*1 -|- 6*2! < 7r/4. Thus, depending on the orientations 
of the inclusions with respect to the separation vector, the interaction between two non- 
circular inclusions can be attractive or repulsive. In general, 2-fold and 4-fold anisotropics are 
nonvanishing and the resulting interaction has more complicated orientational dependence 
as shown in Eq. (2.3). In the case of broken up-down symmetry, the interesting aspect of 
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the interaction in Eq. (2.4), in addition to the orientational dependence of the force, is the 
leading distance-dependence 1/B? rather than The angular dependence in this case 

is of the form cos2(0i + 62) and the interaction is attractive if \9i + < 7r/4 and repulsive 
if 7r/4 < l^^i + ^^2! < 7r/2. In general, transmembrane proteins are asymmetrically embedded 
in the membrane, and they break up-down symmetry of bilayer membranes. Thus, these 
proteins interact with anisotropic 1/E? interaction, which is much stronger at large length- 
scale than the screened electrostatic interaction or the Van der Waals interaction under 
physiological conditions. Consequently, for a distance large compared with a typical protein 
size, the interaction described in this paper will dominate over the electrostatic and the Van 
der Waals interactions. 



Recently, we received the preprint by Golestanian et al. [|T2[ in which they extend the 
calculation of Ref. to the interaction between two rods on membranes. In their work, 
they retain the up-down symmetry for the rods and obtain anisotropic l/R^ interaction 



similar to Eq ( |3.38| ). Also, we'd like to mention the calculation of the short-ranged induced 



interactions between inclusions embedded in fluid membrane by Dan et al. |H| , P^ . They 
find a short-range repulsive interaction with decaying oscillation with period of order the 
membrane thickness when the two halves of a bilayer membrane are allowed to respond sep- 
arately to the membrane inclusions. Consequently, they suggest that in systems where the 
inclusions impose specific contact angles, a meta-stable state with a well-defined separation 
between neighboring inclusions is possible and the minimal energy state of the membrane is 
obtained at a finite inclusion spacing. The competition between this short-range repulsive 
force and the long-range attractive forces discussed by Goulian et. al. and in this paper 
could lead to a preferred separation between membrane proteins. Thus, it may be sug- 
gested that a hydrophilic channel through the bilayer can be formed by a ring of three or 
more transmembrane proteins, which may give some idea about how protein molecules can 
facilitate the passage of ions or molecules into and out of cells. 

We are grateful to Mark Goulian for helpful discussions and for pointing out an error in 
our treatment of the cutoff in an early version of this manuscript. This work was supported 
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FIGURES 




FIG. 1. Proteins make angles 9i measured with respect to the separation vector R. The distance 
between proteins R is taken to be much larger than protein size. 




FIG. 2. Contact angle is measured between the direction of i-th protein and the unit normal 
of the membrane at protein's perimeter. 
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(a) 




FIG. 3. (a) The unit vector m denotes the direction of the protein and N denotes the unit 
normal vector of the membrane at the perimeter of protein (Side view), (b) b is the unit vector 
from the center of the protein to its perimeter (Top view). 




FIG. 4. Non-transmembrane proteins have preferred center-of-mass positions nat at the center 
of the bilayer. ( denotes the position of the protein and h is the membrane height fluctuation field. 
The potential energy has a minimum at the non- vanishing value oi ( — h = ro- 
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